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ABSTRACT 

Consider  the  Weibull  distribution.  The  parameters  of  this  distribution  are  estimated  by  the  maximum  likelihood 
method  and  Bayes  method.  In  this  study  we  present  shrinkage  estimator  between  maximum  likelihood  9 M and  Bayes 

/v 

estimators  0 B that  make  mean  square  error  (MSE)  less  than  other  standard  method.  Using  linear  combination  between 

maximum  likelihood  method  and  Bayes  method  to  obtain  a new  estimator  9 then  simulation  study  will  used  to  compare 

between  shrinkage  estimator.  Maximum  likelihood  estimator  and  Bayes  estimator  to  find  the  best  (less  mean  square  error) 
with  different  sample  size  and  Mathlab  program. 
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1.  INTRODUCTION 


Chiou  and  Han  studied  shrinkage  estimators  of  the  parameter  exponential  distribution;  they  introduced  the  usual 
preliminary  test  estimators  of  the  threshold  parameter  of  the  exponential  distribution  in  censored  samples.  Morris,  Baggerly 
and  Coombes  developed  a new  Bayesian  estimation  procedure  that  quantifies  prior  information  about  two  characteristics, 
yielding  a nonlinear  shrinkage  estimator  with  efficiency  advantages  over  the  MLE.  In  1996  Singh  and  Raghuvanshi  studied 
the  problem  of  estimation  of  variance  in  exponential  density  when  a prior  point  estimate  is  available.  In  2008  Baklizi  and 
Ahmed  studied  three  classes  of  point  estimators,  namely,  the  unrestricted  estimator,  the  shrinkage  estimator  and  shrinkage 
preliminary  test  estimator  and  mean  squared  errors  are  derived  and  compared. 


In  this  paper  we  studied  shrinkage  estimator,  0 between  Bayes  and  Maximum  likelihood  estimators.  Using 

simulation  study,  we  compared  between  estimators  shrinkage  estimator.  Maximum  likelihood  estimator  and  Bayes 
estimator  of  parameter  Weibull  distribution  to  find  the  best  depend  on  less  mean  square  error  MSE. 

The  estimators  of  parameter  Weibull  distribution  was  estimated  (see  A1  Omari,  Alkutubi  and  Akma2010), 
supposed  be  a random  sample  of  size  n with  distribution  function  and  probability  density  function.  Using 

Bayesian  and  Maximum  likelihood  (MLE)  methods  to  obtain  these  estimators  The  probability  density  function  of 


Weibullcase  is  f ( t ; 9 , p ) 


^tp~l  exp( 

9 


t p 

9 


) , and 


On 


X 

_i=i is  the  Bayes  estimator.  The  mean 

n — 1 


square  error  of  Bayes  estimator  is  given  by  expected  value  and  variance. 


such  that  MSe(gb ) = ” + 1 G 2 
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E {eR)=  —"—0  Var(dB)=  n e- 
V s ’ n - 1 In  -1) 

n 

Also  the  maximum  likelihood  estimator  of  exponential  distribution  is  given  by q = <- , And 

n 

0 2 

MSE  (6  ) = — , is  the  mean  square  error  of  maximum  likelihood  estimator. 

n 

2.  MATERIALS  AND  METHODS 

2.1  Shrinkage  Estimator  §]  between  Bayes  and  MLE 

In  this  section,  we  obtain  shrinkage  estimator  6 j Weibull  distribution  from  Bayesian  and  maximum  likelihood 
estimators  by  linear  such  that, , 

#1  = PAm  + (!  - Pl)0B 

The  value  of  Pi  which  minimizes  MSE(6{)  is,  p { = — — + 2 n + 

4 n 2 - n + 2 


Q 

Then  the  shrinkage  estimator  1 between  Bayes  and  MLE  of  Weibull  distribution  is  given  by 


= ( 


n " + 2 n + 1 
4 n 2 — n + 2 


)0M  + (1-  ( 


n +2/7  + 1 
4 77  2 - 77  + 2 


))0j 


2.2  Shrinkage  Estimator  § Between  Q and  MLE 


We  can  obtain  the  shrinkage  estimator  0 , for  Weibull  distribution  from  0 x and  MLE  depending  on  linear 
combination  to  get  the  following  equation, 

^21  — Plfi\  + (1  — 


To  find  the  value  of  Ell  which  minimizes  MSE  ($2i  ) > we  follows 

Q2\-Q= iPiA  + (i  - Pu)6m  ] ~ e 

Taking  expected  valued, 

E(02l  -ef  = p^E^  -Of  +(1-  p2l)2E0M  -Gf  + plx{\-plx)E(Px-G\eu-G\ 

Then 

MSEW2])  = p22\MSE{6x ) + (1  - p2l)2 MSE(0m ) + 2p2l(l-  p21)[E(0pM ) - E{9X0)  - E0Md)  -E02], 


Impact  Factor  (JCC):  2.6305 
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Let  d MSE { 9 ^ , implies  the  value  of  which  minimizes  MSE  (9  ^ ) is 


dP 


Pn  = 


MSE  (pM  )-  E (epM  )+  e{bp  )+  E (eM  e )-  e q?): 


E ( 9 J 


MSE  (9l  )+  MSE  {§ M )-  2 E (<9,<9M  )+  2 E (di9  )+  2 E (S M 9 )-  2 E (fl  2 ) 
4 zz  6 - 6 zz  5 + zz  4 + 2 zz  3 - 6 zz  2 -In 


, where 


v zz  5 - 10  n 5 + 12  zz  4 - 18  n 3 + 9 n 2 - 3 zz y 


E(oM)  = e 


MSE  (9X)  = 


8 /z 8 — 2 zz 7 — n 6 + n 5 — 6 n 4 — 4 zz 3 + 6 zz  2 + 2 n 
2 n 8 - 29  zz 7 + 25  zz 6 - 24  n 5 + 33  n 4 - 10  zz 3 + 16  n 2 - 4 n + 2 


and 


MS£  (0M)  = — P 21  ~ 

n , Then 


X 


21 


^2 


1 , where 


6 zz 7 - 14  zz 6 + 19  zz 5 - 12  zz 4 + 16  zz  ' - 12  zz  ’ + 2 zz 


y,  L (n  - 2 n 2 + 10  zz  3 - 12  zz  4 + 13  zz  5 - 12  zz  6 - zz  7 

14  zz 8 - 22  zz 7 + 27  n 6 - 26  zz 5 + 20  zz  4 - 2 zz 3 + zz 2 

V / 


6> 


and. 


_y  2I  = [(57  ft15  + 29  n 14  + 4540  n 13  - 361  n 12  - 670  n 11  + 2231  n 10 

- 172  n 9 + 142  n 8 - 565  n 7 - 66  n 6 + 19  n 5 - 211  n 4 + 121  n 3 

- 40  « 2 + 21  ft  - 4)5-  (2ft  — 10  ft  2 — 42  ft  3 - 130  n 4 + 322  ft  5 

- 652  n 6 + 1342  zz  7 - 2431  n 8 + 2784  zz  9 - 2349  zz  10  + 2958  zz  11 

- 1226  zz  12  + 306  zz  13  - 135  zz  14  +27  zz  15  )]  <9  2 + 


r z 

(-  2 zz  7 — 6 zz  6 + 

31  zz  5 

- 29 

zz  4 + 

33  zz  3 - 4 n 2 - 3 zz 

)] 

5 zz  8 - 17  zz  7 + 13 

n 6 - 

14  zz 

5 + 12 

zz  4 - 3 zz  3 + zz  2 

Then  the  shrinkage  estimator  of  Weibull  distribution  02\ 


between 


*1 


and  MLE  is  given  by 


elx  = + (i  - 

y21  y21 

^ fj 

2.3  Shrinkage  Estimator  22  Between  ! and  Bayes 

In  the  same  way  in  above,  we  can  get  the  shrinkage  estimator  ® 21  of  Weibull  distribution  between  0 x and  Bayes 

such  that 


9 22  P 22^  1 “E  ( 1 P 22  ")9g 
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The  value  of  P22  which  minimizes  MSE  (022 ) as  follows 

MSE  (0B)~  Eje ig)  + E(8B8)  - E{92) 

Pn  MSE  (0,)  + MSE  (0B)~  2E(010b)  + 2 E(dB0)  + 2 E(dx0)  - 2E{62) 

Depending  on  MSE{0B ) , MSE ( 0i ) , E(0B)  and  E(0^  ) we  can  obtain  P 22 , such  that 


P22  , Where  A22 , )?22 , W22  and  Z22  is  given  by 

^22  ^22 


X 

22 

60  n 16 

- 570  n 15  + 2235 

n 14  - 5642 

13 

n 

+ 9721 

12 

n 

52310 

11 

n 

+ 16940 

n 10  - 18421  n 9 + 

16321 

n 8 - 

98542  n 7 

+ 

56542 

n 6 - 

39642  n 

5 + 17890  n 4 - 553  n 3 + 167 

n 2 — 

12 

n + 6 

y 

22  = 

IO  16 

12  n 

- 123  n 15  + 743  n 14 

- 2656 

n 13 

+ 

4643  n 

12 

6750 

n 11 

+ 9378 

n 10  - 13674  n 9 

+ 9321 

n 

8 

7033 

n 1 

+ 4476 

n 6 

- 2287 

n 5 + 926  n 4 - 390 

n 3 + 

77  n 

2 _ 

9 n + 4 

w 

22 

= 130  n 25  - 

- 1787 

77  24  + 

11652 

77  23  - 58750  77  22  + 

148702  77  21 

- 338904 

n 20  + 599433 

77  19  - 

587526 

77  18  + 687532 

17 

77  - 

1223842 

n 16  + 534329 

n 15  — 

787456 

77  14  + 590818 

77  13  — 

619083 

n 12  + 523976 

77  11  - 

390395 

77  10  + 254273 

77  9 - 

213304 

n 8 + 56064 

« 7 

- 25698 

77  6 + 

106114  77  5 - 

6591 

n 

4 + 1775 

n 3 

- 179 

77  2 + 

34  77  - 2 

z 

22 

= 18  n 24 

- 169 

77  23  + 

944  77 

22  - 3238  77  21  + 

7754  77  20 

- 15453 

n 19  + 27895 

« 18 

- 28800 

77  17  + 

54302  77  16  - 45678 

77 15 

+ 57898 

n 14  - 56509 

77  13 

+ 45748 

12 

77  - 

45433  77  11  + 34592 

10 

77 

- 37894 

n 9 + 24551 

n 8 - 

- 15138 

77  7 + 8060  77  6 - 3622  n 5 

+ 1336  n 

4 

- 387  n 3 + 79 

n 2 — 

9 77  + 1 

Then  the  Shrinkage  Estimator  0 Between  0 and  Bayes  is  given  by 

822  = (^-x  M^-)8l  + (1  - ( ^22  x VV22  ))  0B 

y 22  z 22  y 22  z 22 

3 SIMULATION  STUDY 

In  this  study,  we  chooses  samples  sizes  n=30,  60,  90,  with  parameter  value  6 = 0.5,  1,  1.5  and  R=1000  of 
replication.  Using  mean  square  error  (MSE)  to  compare  between  all  estimators  6y , 6^x  , 021 ,0sand0M  , where 

t (a,  - « J 

MSE  (8  ) = , 

R 

The  simulation  program  is  written  by  Matlab  program.  And  the  results  are  introduced  and  tabulated  in  Table  1 for 


Impact  Factor  (JCC):  2.6305 


NAAS  Rating  3.19 
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the  MSE  of  all  estimators  for  all  sample  sizes  and  0 values  respectively. 


Table  1:  MSE  of  Estimator  Exponential  Distribution 


N 

e 

e 21 

°22 

@ B 

0.5 

0.0324 

0.0320 

0.0321 

0.0312 

0.0312 

30 

1 

0.0322 

0.0318 

0.0319 

0.0321 

0.0325 

1.5 

0.0310 

0.0307 

0.0309 

0.0309 

0.0311 

0.5 

0.0226 

0.0220 

0.0221 

0.0222 

0.0222 

60 

1 

0.0218 

0.0217 

0.0218 

0.0219 

0.0219 

1.5 

0.0217 

0.0216 

0.0217 

0.0220 

0.0221 

0.5 

0.0134 

0.0126 

0.0129 

0.0131 

0.0135 

90 

1 

0.0088 

0.0085 

0.0086 

0.0087 

0.0089 

1.5 

0.0050 

0.0047 

0.0048 

0.0048 

0.0049 

From  this  Table,  the  shrinkage  estimator  0 2\  is  the  best  shrinkage  estimator  from  all  estimators  and  for  all 
sample  size  and  parameter  value. 


4 CONCLUSIONS 


The  generator  shrinkage  estimator  0 2\  is  the  best  estimator  (less  mean  square  error).  The  effects  of  sample  size 
on  the  mean  square  error  of  all  estimators  refer  to  the  MSE  decreases  as  n increases. 
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